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Abstract 

We analyze an adaptive finite element/boundary element procedure for 
scalar elastoplastic interface problems involving friction, where a non- 
linear uniformly monotone operator such as the p-Laplacian is coupled 
to the linear Laplace equation on the exterior domain. The problem 
is reduced to a boundary/domain variational inequality, a discretized 
saddle point formulation of which is then solved using the Uzawa al- 
gorithm and adaptive mesh refinements based on a gradient recovery 
scheme. The Galerkin approximations are shown to converge to the 
unique solution of the variational problem in a suitable product of L p - 
and L 2 -Sobolev spaces. 

1 Introduction 

Consider the following transmission problem on a bounded Lipschitz domain 
Q C R n : 

-div(e(|Viii|)Vtii) = / inft, -Au 2 =0 in Q c , 
g(\Vui\)d u ui - d v u 2 = t on d£l,ui - u 2 =u on T t , 
(1) -^(iVniD^u^uo + u 2 - ui) + g\(u Q + u 2 - ui)| =0, 

|e(|Vui|)<?„ui| < g on r s . 

. . f o + o(l) ,n = 2 
n2(x) = | 0(|x| 2 -),n>2 • 

Here g(t) denotes a function g(x,t) G C(0 x (0, oo)) satisfying 

0<g(t)<g*[t\l + t) 1 - s ] p - 2 , 
\g(t)t - g(s)s\ < g*[(t + s)\l + t + s) 1 - 5 ]^^ - s\ 
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and 

g(t)t - g(s)s > g*[(t + s) 5 (l + t + s) l - & f- 2 {t - s) 

for all t > s > uniformly in x G f2 (5 G [0, 1], g*, g* > 0). The interface 
d£l = r s U T t is divided into the disjoint components T s and T t ^ 0, and 
the data belong to the following spaces: 

/ G n G PH' 2 (dtt), t G wM' 2 (dft), g G L°°(r s ), a G R. 

As usual, the normal derivatives are understood in terms of a Green's for- 
mula, and it is convenient to set a = for n > 2. In two dimensions one 
further condition is required to enforce uniqueness: 

(2) / /+<to,l> = 0. 

Jn 

We are looking for weak solutions (ui, U2) G W 1,P (Q) x W^(Q C ) whenp > 2. 

A typical example is given by g(t) = [t s (l + t) 1 ^ 5 ]^ 2 , 5 G [0, 1], with the 

p-Laplacian corresponding to the maximally degenerate case 5=1. 

In this article we use layer potentials for the Laplace equation on f2 c to 

reduce the system to a uniquely solvable variational problem on W 1,P {^1) x 
I 2 

Wq' (r s ). The main idea of our theoretical analysis is simple: Because 
the traces of W^ 1,p (f2)-functions are continuously embedded into Wz' 2 {dVL) 
for p > 2, the quadratic form (Su, u) associated to the Steklov-Poincare 
operator is accessible to Hilbert space methods whenever it is defined. In 
this slightly weaker setting, Friedrichs' inequality (Prop. [T]) allows to re- 
cover control over the LP -norms in the interior, and as a consequence the 
full variational functional associated to the above equations is coercive in 
W l ' p (n). 

In the numerical part we present a model problem, which shows singu- 
larities resulting from the given boundary data, as well as from the change 
of boundary conditions, leading to a suboptimal convergence rate for uni- 
form mesh refinements. We also present a Uzawa solver to deal with the 
variational inequality. 

With the help of a Korn inequality (Prop. [2]), our method easily carries 
over to transmission problems in nonlinear elasticity, e.g. Hencky materi- 
als in f2 coupled to the Lame equation in Q c . A generalization to certain 
nonconvex energy functionals will be discussed elsewhere [Tj- 

The outline of the article is as follows: Section[2]recalls some properties of 
L p -Sobolev spaces and introduce a family of quasinorms adapted to the con- 
sidered class of operators. In the following section [3] we introduce the bound- 
ary integral operators and derive our variational formulation. Section H] is 
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dedicated to the existence and uniqueness of our model problem. The dis- 
cretization of our problem is derived in section [5l as well as the a-priori error 
estimates. In section [6] our a-posteriori error estimator is presented and its 
reliability proven. Finally, in section [7] we present the Uzawa-solver and two 
numerical examples, clearly underlining our theoretical results. 

2 Preliminaries 

Let £1 be an open subset of M n with Lipschitz boundary d£l. Set p' = 
whenever p E (1, oo). 

Definition 1. The Sobolev spaces W^ P (Q), k E No, are the completion of 
C^(J7) with respect to the norm \\u\\ W k, P ^ = \\u\\k :P = \\u\\p+Y^\ y \=k W^uWp. 
The second term in the norm will be denoted by \u\yyi, P m\ = Let 

W^ Kv '{Vt) = (W k >P(n))' and W~ k 'P\n) = (w o fc,p (n)Y. W^p^idn) de- 
notes the space of traces of W l,p ($V) -functions on the boundary. It coincides 

with the Besov space Bp. p p (dQ) as obtained by real interpolation of Sobolev 
spaces J77] /. and one may define W s ' p (dQ) = Bp, p (dft) for s E (—1, 1). 

Remark 1. We are going to need the following properties for bounded dVl 

EZJ-- 

a) All the above spaces are reflexive and (W s ' p (dQ))' = W s,p (dfl). 

b) For p = 2 they coincide with the Sobolev spaces H s . 

c) W X ~ ^ W^' 2 (dQ) for p>2. 

d) If dO, is smooth, pseudodifferential operators of order m with symbol in 
the Hormander class S™ (<9r2) map W s,p (d£l) continuously to W s ~ m ' p {dVt). 
For Lipschitz dQ, at least the first-order Steklov-F 'oincare operator S of the 
Laplacian on Q c is continuous between Wz ,2 (dVl) and W~z ,2 (dQ) J^j. 

e) Points a) to d) imply that the quadratic form (Su,u) associated to S is 

well-defined on W l ~p' p {d£l) if p > 2. S being elliptic, the form cannot be 
defined for p < 2 even if d£l is smooth. 

Uniform monotony will be shown using a variant of Friedrichs' inequality. 

Proposition 1. Assume f2 is bounded and that T C <9f2 has positive (n— 1)- 
dimensional measure. Then there is a C > such that 

\\u\\ P < C(\\Vu\\ p + Hr|Ui(r)) for all u E W lj, (Q). 



3 



Proof. We apply an interpolation argument to the well-known Friedrichs' 
inequality 

1 



\u - un\\p < C||Vu||p, u n = j— j u, 



n 

on W lj, (Q) (see e.g. [10]). Let L : W^fi) -» L p (ft) be the rank-1 operator 
Lu = / r ujr and 7 the inclusion of W 1,p (f2) into L p (f2). Then I — L : 
W 1 * (O) -> L p (ft) is bounded and 

||« - Lu|| p = - L)(u - uq)\\ p < \\I - L\\\\u - un||i,p < C||Vit|| p 

for all u € VF 1,p (f2). The assertion follows. □ 

Let u(x,y) = (\x\ + \y\) s {l + \x\ + < S < 1. In addition to the 

above norms, the following family of quasi-norms will prove useful: 

Definition 2. For v,w € VK 1,P (S7) and G No, define 

i 

\v\ {k>w , P ) = I J u{Vw,D k vf- 2 \D k v\^ 2 
where \D k v\ 2 = S| 7 |=fc l«9 7 v| 2 . 

Remark 2. a) If p > 2, the (l,w,p)-quasi-norm can be estimated from 
above and below by suitable powers of the W 1,p -seminorm Jfjjj: 

Mi, P < \ v \ 2 (i,w,p) < c(Mi,p> \ w \i,p)\ v \i,p- 

b) In the nondegenerate case 5 = 0, we have \v\\ 2 < M(i W pV 

c) The following inequality is useful for computations with quasi-norms: 

Xp < maxjr 1 , e^-rijia?- 1 + A) P '~ 2 A 2 + e(a + p) v ' 2 p 2 
for A, p, a > and e > 0. 

The results of this paper easily generalize to the systems of equations 
describing certain inelastic materials. In this case, Lemma [1] has to be 
replaced by the following Korn inequality: 

Proposition 2. Assume C W 1 is a bounded Lipschitz domain and V C dQ, 
has positive (n — 1) -dimensional measure. Then there is a C > such that 

Nil* < C(\\e(u)\\ p + |Mr|| L1( r)) for all u E (W^(n)) n . 
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Proof. The IP -version ||n||i iP < C(||e(n)||p + \\u\\ p ) of Korn's inequality is 
well-known (see e.g. [5]). Assume the assertion was false. Then ||e(tt n )||p + 
ll u n|r||Li(r) — \ f° r some sequence in W 1,p (0,) normalized to ||n n ||i jP = 1. 
By the compactness of W l,p (Vl) L p (£2), we may assume u n to converge in 
L p {0). The cited variant of Korn's inequality shows that u n is even Cauchy 
in W ' p (il), hence converges to some no with ||e(no)|| p = ||^o|r||Li(r) = 0- 
The kernel of e consists of skew-symmetric affine transformations Ax + b, 
A = —A T . As dimker A = n mod 2, no cannot vanish on all of the (n — 1— 
dimensional) V unless uq = 0. Contradiction to ||tto||i, p = 1. □ 



3 Variational Formulation and Reduction to dfl 

We continue to use the notation from the Introduction and mainly follow 
[9J- Fix some p > 2 and, for q(t) = J* * sg(s) ds, let G{u) = f n q(\Vu\) with 
derivative 

DG(u,v) = (G'u,v) = [ g(\Vu\)VuVv (u,v G W^ p (n)) 

Jn 

and j(v) = J* r g\v\, v G L (F s ). G is known to be strictly convex and G 1 : 

— > (W 1,P (Q)) bounded and uniformly monotone, hence coercive, 
with respect to the seminorm | ■ There is some «g > such that for all 
u,v & W l,p {Vt) 

{G'u — G'v,u — v) > ac\u — v\\ and lim — : — p — - = oo. 

The naive variational formulation of the transmission problem ([1]) mini- 
mizes the functional 

$(«i,n 2 ) = G(ux) + - I |Vn 2 | 2 - / fui- (t ,u 2 \dn)+j({u2-ui + uo)\r s ) 
1 Jq c Jn 

over a suitable convex set. 

Lemma 1. Minimizing over the nonempty, closed and convex subset 
C = {(u!,u 2 ) G W^(n) x W£ C 2 (Q C ) : (m - u 2 )\r t = n , n 2 G C 2 }, 

£-2 = {v G Wl^(VL c ) : Av = 0in W~ 1,2 (Q C ) + radiation condition at oo}, 

is equivalent to the system |7]j in the sense of distributions if g G C l (Q. x 
(0,oo)). 
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Proof. C is apparently convex. A similar argument as in Remarks 2 and 4 of 
PQ shows that C is closed and nonempty. The proof there almost exclusively 
involves the exterior problem in C 2 and only requires basic measure theoretic 
properties of W 1,2 (Q), which also hold for W 1,P (Q). Finally, repeat the 
computations of [S] to obtain equivalence with ([!]). □ 

To reduce the exterior problem to the boundary, we are going to need 
the layer potentials 

V<p(x) = / 4>( x ') log\x — x'\ dx', 



TT 



an 



JC(ft(x) = / 4>{x') d u , log \x — x'\ dx' , 

Jan 

K,'4>{x) = / 4>{x') dy x log \x — x'\ dx', 

n Jan 

W4>(x) = — d u I Mx') d u , log \x — x'\ dx' 
t Jan 

associated to the Laplace equation on Q c . They extend from C°°(<9f2) to a 

yy ^1 J on the Sobolev space W^' 2 (dQ) x W~z' 2 (dfl). If 

the capacity of dSl is less than 1, which can always be achieved by scaling, V 
and W considered as operators on W~z ,2 (dVL) are selfadjoint, V is positive 
and W non-negative. Similarly, the Steklov-Poincare operator 

s = w + (1 - x:')v -1 (i - k) ■ w^ 2 (on) c w-^ 2 (dn) -> w^' 2 (an) 

defines a positive and selfadjoint operator (pseudodifferential of order 1, if 
d£l is smooth) with the main property 

d„U2\dn = -S(u 2 \aQ - a) 

for solutions u 2 £ C 2 of the Laplace equation on J7 C . By Remark [I] e), S 
gives rise to a coercive and symmetric bilinear form (Su,u) on W^' 2 (dVL) 
and, in particular, a pairing on the traces of W 1,p (ft) if and only if p > 2. 

Using the weak definition of d u \an, S reduces the integral over Q c in $ 
to the boundary: 

/ |Vn 2 | 2 = -(d u u 2 \dn,U2\an) = {S(u 2 \an - a),u 2 \ d n) for u 2 E C 2 . 
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Easy manipulations allow to substitute U2 by a function v on T s (cf. [9]): 
Let 

jyi 2 (r s ) = {u g W2. 2 (5si) : supp u c r s }, x p = w^in) x W5' 2 (r s ) 

and (u,v) = (ui—c,UQ+U2\dn — ui\dn) G X p for a suitable c£R. Collecting 
the data-dependent terms in 

X(u,v) = {t + Su ,u\ 9n + v) + / /u 

leads to 

$(«i, w 2 ) = G(u)+-{S(u\ dn +v), u\ dn +v)-X(u, v)+j(v)+~(Su , u )+(t , u ) 

The first three terms on the right hand side will be called J(u,v). 

Lemma 2. Minimizing $ over C is equivalent to minimizing J + j over 
the nonempty closed convex set D = {(u, v) G X p : (S(u\ofi + v — uo), 1} = 
ifn = 2} 

Proof. As in [9]. The main additional observation here is that the substitu- 
tion v = uq + U2\oq — ui\dQ indeed defines an element of W^' 2 {T S ), because 
u ,u 2 \dn G W^' 2 (dn), Mi| 9 n G W l '^ p {dVt) C W^' 2 (<9ft) by Remark^and 
u|r t = 0, if (ui,U2) G C. □ 

4 Existence and Uniqueness 

Minimization of J + j over D translates into the following variational in- 
equality: Find (u, v) G X p such that 

(G'u, u - u) + (S(u\ dn + v),(u-u)\qq + v- v) +j(v) - j(v) > X(u -u,v-v) 
for all (it, v) G X p . Note that D has been replaced by X p . 

We now prove the crucial monotony estimate: 

Lemma 3. The operator in the variational inequality is uniformly monotone 
on X p . There exists an a = a(C) > such that for all \\u,v\\x, \\u, v\\x < C 

a (h-u\\w*{n) + \\ v -*\\~i,3pJ ^ (G'u-G'u,u-u) 

+ (S((u -u)\ dn + v - v), (u -u)\an + v- v). 
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Proof. Recall the monotony estimate for G' from Section [3j 

(G'u — G'u, u — u) > ac\u — u\\ . 
The triangle inequality and convexity of x p imply 

lit) — < (ll(tt — tt)|r„ + v — v II i, + \\(u — tt)|rj| 1, ) p 

< 2^(ii( fl -„)| r , + i -„ii;, l(r _ ) + ii( 4 -„)| r ,ii;,^ )) . 

Using VF 1 p' p (r s ) Ty2' 2 (r s ) as well as the boundedness of the trace 
operator, 

2 1_p \\v - f||Li - /3 \\u-u\\^ ltPm <\\(u-u)\ Ts +v -vf x 

follows for some f3 > 1. Let 

K = {(u,v,u,v)eXPxXP:\\(u-u)\ dQ + v-v\\ wh2{an) <2pC} 

and < e < [3~ l . Since S 1 is positive definite on W^' 2 (dQ), we ob- 
tain from Priedrichs' inequality for (u, v, u,v) G K or, in particular, if 

||«,u||x, \\u,v\\x < C: 

{G'u - G'u,u - u) + (S((u - u)\dn + v — v), (u - u)\dn + v - v) 
>\ u -u\l p+ \\( u - u )\ m + v- v f wh2 ^ 

Z\u-<, P + \\(u-u)\ 9 n + v-v\\^ h2 ^ 

>\u- u\l p + e\\(u-u)\ Ts +v- + \\(u - tOlrJ^^ 

> ||« - u\\ p wl , p{n) + e\\(u-u)\ r .+v- w|l^, 3(I , a) 

> (1 - s(3) \\u - uf whpm + 2^ p e \\v - v\\^ h2{Ts . 

Uniform monotony on all of X p is shown similarly, but on the unbounded 
complement (X p x X p ) \ K the exponents p on the left hand side have to 
be replaced by 2. □ 

Theorem 1. The variational inequality is equivalent to the transmission 
problem |ip and has a unique solution. 

Proof. We repeat the computations in [9] to get the equivalence with the 
minimization of J + j over D, and hence with ([1]). Existence and uniqueness 
follow from Lemma El e.g. by applying [12], Proposition 32.36. □ 
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5 Discretization and Error Analysis 

In order to avoid using S = W + (1 — X7)V _1 (1 — fC) explicitly, the numerical 
implementation involves a variant of the variational inequality 

(G'u, u - u) + (S(u\qq + v) , (u - u) \an + v - v) +j(v) - j(v) > X(u -u,v-v) 

in terms of the layer potentials. Our a posteriori analysis is therefore based 
on the following equivalent problem: Find (u,v,(f>) G X p x W~i ,2 (d£l) =: 
Y p , such that 

{G'u, u-u} + {W(u\ 9 q + v) + (K! - 1)0, (u -u)\aa + v- v) 

+ j(v)-j(v) > (t + yVu ,(u-u)\ dn + v-v)+ / f(u-u), 

Jn 

(0, Vj> + (1 - K)(u\aa + v)) = (0, (1 - K.)uq) 
for all (u, v, <j)) G More concisely, 

B(u, v, 4>; u — u, v — v, 4> — 4>) + j(v) — j(v) > A(u — u, v — v, <p — <j>) 
with 

B(u,v,4>;u,v,4>) = (G'u,u) + (W(u\ 9n + v) + (K! - l)<f>,u\gn + v) 

+(0,V0+(l-XOHai + «)), 
A(u,v,^>) = (t + Wu ,u\ dn +v) + / /u+ (0, (1 -TC)uq). 



The more detailed a priori and a posteriori error analysis requires a few 
basic properties of the quasi-norms [6]. 

Remark 3. a) The continuity and coercivity estimates can be sharpened: 
For allu,v G W l *{Q,) 

(G'u - G'v,u- v)<\u- v\f 1>Ujp) < (G'u - G'v,u- v). 

b) There is 6 > such that for all e G (0, oo) and all u,v,w G W 1,p (£i) 

| (G'u - G'v,w}\ < e\u - v\ 2 (1 ^ p) + e~ e \w\ 2 {l ^ p) . 

Lemma 4. For all (u,v,(f)),(u,v,<fi) G Y p we have 

\u - u\ 2 x ~ } + \\(u - u)\dn + v - v\\ 2 t + \\r]\\ 2 ! 

^ \u - u\% ,, n) + \\(u - u)\ dn + v - v\\ 2 i „ +U-6\\ 2 lo 

< B(u, v, eft; u — u, v — v, rj) — B(u, v, <fi; u — u, v — v, if), 

where 2rj = <p — <P + V _1 (l — K)((u — u)\qq + v — v). 
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Proof. The right hand side of the identity 

B(u, v, 4>; u — u, v — v , rj) — B(u, v, 4>; u — u, v — v , rj) 
= (G'u - G'u,u- u) + \(W((u - u)\on + v - v), (u - u)\on + v - v)) 
+ - u)|«i + v - v), (u -u)\ dn + v- v)) + \(V{4> -4>)J- <f>). 

is, up to a constant, larger than ||u — u, v — v, ij> — 0||^ YP y Furthermore, 

□ 

Let {Th}h & i a regular triangulation of O into disjoint open regular trian- 
gles K, so that = Ukgt^ Each element has at most one edge on dCl, 
and the closures of any two of them share at most a single vertex or edge. 
Let hx denote the diameter of K € 7^ and px the diameter of the largest 
inscribed ball. We assume that 1 < max^ e r h t£ < R independent of h and 
that h = m&XK£T h hx- £h is g° m g to be the set of all edges of the triangles 
in Th, D the set of nodes. Associated to 7/j is the space 
of functions whose restrictions to any K € 7^ are linear. 

i 2 

50 is triangulated by {Z € £^ : Z C d£i}. W h 2 ' (dQ) denotes the corre- 

i,2 

sponding space of piecewise linear functions, and W^ 2 ' (r s ) the subspace of 

_l 2 i 

those supported on T s . Finally, VF fc 2 ' (dQ) C TV" 2 > 2 (dQ). 



We denote by i h : W^Sl) W^ft), j h : W ft * ,2 (r<) ^ W^ 2 ^, 



_1 2 x 

and fe/j : W. 2 ' (dQ) <^-> 1F~2> 2 (<90) the canonical inclusion maps. Set 
X£ = W^' p (0) x ^' 2 (r s ), We denote by i h : I^' p (0) ^(fl), ^ : 
^ 2 ' 2 (r s ) ^i 2 (r s ) and : W/^' 2 (dO) VF"5- 2 (aO) the canonical 
inclusion maps. Set X v h = W^ p (n) x w\ ,2 (T S ), 

Sh = \{W + (/ - K')k h {kiVk h T x kt,{I - K)) 



and 



K{uh,v h ) = (t + S h Uo,u\gQ + v) + / /u/j. 



As is well-known, there exists Zio > such that the approximate Steklov- 
Poincare operator Sh is coercive uniformly in h < ho, i.e. (ShUh,Uh) > 
a s\Wh\\ 2 w i 2 ^ with as independent of h. 
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The discretized variational inequality reads as follows: Find (uh, Vh, 4>h) 6 
Yfi such that 

B(u h ,v h , <f> h ; u h -u h , v h -v h , 4>h-(f)h)+j(vh)-j(vh) > A(u h -u h , v h -v h , 4>h-4>h) 

for all (uh,Vh, <fih) £ Y?. Repeating the arguments from the previous section, 
one obtains a unique solution to the discretized variational inequality. 

Theorem 2. Let (u, v, (p) € Y v , v^, (frh) G be the solutions of the con- 
tinuous resp. discretized variational problem. The following a priori bound 
for the error holds uniformly in h < h$: 



\\u -U h ,V- V h ,<t>- <t>h\\YP 

< [« - + IK* - + v ~ v h \\ 2 w i, 2{m) + M- kf w ^ {m) 

< inf \\u-u h ,v-v h ,$-<t> h \\ Y p + \\v-v h \\ L 2( Ts y 

(u h ,v h ,<f> h )£Y£ 

Proof. Let (u,v,cp) € Y p , (u^, Vh, (frh) € Y£. Lemma H] and the variational 
inequality imply 

I* " fifclfx,*,) + m - fi*)|an + 6 - «fcll^^ (80) + II* " ^^-^n) 

< jB(tt, v,(f>;u- u h , v-v h ,4>- cj) h ) - B(u h , v h , (f> h ; u - u h , v - v h , <f> - (f> h ) 

< B(u, v, 4>; u, v, (f>) - A(u - u, v - v, <f> - $) + j(v) - j(v) 

+ B(u h ,v h , 4> h ; u h , v h , <t>h) ~ Huh ~ u h ,v h - v h , (p h - <j) h ) + j{vh) - j{vh) 
- B(u h ,v h , (f> h ; u, v, 4>) - B(u, v, cj>; u h ,v h , 4> h ) 

Setting (u,v,cj)) = [uh,Vh,4>h) an d adding 0, the right hand side turns into 

B(u, v, 4>; u h u, Vh v, 4> h - 4>) - Mu h -u,v h - v, 4> h - 4>) + j(v h ) - j(v) 
+ B(u, v,cj);u- u h ,v -v h ,(f)- <j> h ) - B(u h , v h , (f) h ; u - u h , v - v h , 4> - <f>h)- 

We first consider the friction terms: 



j(vh)-j(v)= g(\v h \-\v\)< / g(\v h -i>\) <\\g\\ L 2^ a) \\v h -v\\ L 2 {rs) . 

The last two terms are bounded using Remark [3b and Cauchy-Schwarz: 

(G'u-G'u h ,u-u h ) < e\u - Uhlfi^p) + - Uh\fi >ii)P ), 

< e\u h - u\\x Ap ) + s~ e C(\u\ 1:P , \u h \i, p )\u h - u\j iP 
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for sufficiently small e > 0. We may replace C(|tt|i )P , |w/i|i, p ) by an honest 
constant noting that the coercivity of our functional gives an a priori bound 
on an d that we can restrict to those uh satisfying ||ttfc||wi>i>(«) — 

2||«||^i, P (n). Moreover, 

(W((u - u h )\an +v-v h ) + (l- JC')(4> - 4> h ), (u - u h )\an + v -v h ) 
< ell (u — Uh Uo + v — Vh II 2 i„ + e\\d> — 6 h II 2 x „ 

+ e^ 1 \\u — Uh)\\ 2 i, + e" 1 llv — Vh\\ 2 i„ , 

and 

(<£ - <t>h,V{4> - 4>h) + (l - - "jOlan + u - %)) 



2 



1 



Substituting (u,v,(f>) = (u/i,{),0) and (u,v,(f>) = (2u — Uh, v, 0) into the 
variational inequality on Y p and using that also the <fi part is really an 
equality, the remaining two terms reduce to 

(-to - Wu + W(u\ an + v) + (£' - 1)0, v h - v) 

= -(t - S(u\ 9 n +v - u ), v h - v) 

= -(g(\Vu\)d u u,v h - v) < ||5llL 2 (r s )ll^ - *llL 2 (r a )- 

Applying these various estimates to the terms of the right hand side, the 
assertion follows from 

\\u-Uh,v-Vh,$-$h\\Yv % l^"^l?i,«,P) + ll^"^)l^ + *~^ll^,2 (an) + l^~^ll^-^2 (9n) 
as in Lemma El □ 

In the nondegenerate case 5 = 0, we essentially recover the estimates for 
uniformly elliptic operators from [U [9] . 

Corollary 1. For 5 = 0, we obtain 

\\u-u h ,v-v h ,(f>-(t> h \\ Y 2 < inf \\u-Uh,v-v h ,$-<p h \\ YP +\\v-v h \\ L 2(r a ) 

{u h ,v h ,<t>h)&Y£ 

uniformly in h < ho 

Proof. Use (2)3) to estimate |% — u\nu,p) in Theorem [2] from below. □ 
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6 A posteriori error estimate 

Denote by 

(e, e, e) = (u - u h ,v -v h ,(f>- $ h ) € Y p 

the error of the Galerkin approximation, and let 2u = e+V _1 (l— ^C)(e|an+e). 
Our basic a posteriori estimate is the following. 

Lemma 5. For all (e^, %, Uh) £ Yfi 

\^(i,u, P ) + Iklflo + zf w ^ m + Ml^-i. a(80) 

< A(e - e h , e - e h , u - u h ) + j(e h + *5 ft ) - 
- B(u h , v h , (j) h ; e - e h , e - e h , u - u h ) 

f(e - e h ) - {G'u h , e - e h ) + / g(\e h + v h \ - \e + v h \) 



- (u - u h , V(f> h + (1 - fC)(uh\dn + v h - u )) 

+ (t - W(uh\an + v h - u Q ) - (K 1 - l)4> h , (e - e h )\an + e-e h ). 

Proof. Lemma HI the continuous and the discretized variational inequality 
imply 

\u - u\n f, n) + - u)\dn + v - v\\ 2 i 9 +||^-0|| 2 i„ 

< J3(«, v,4>;u- u h ,v - v h , v) - B(u h , v h , (f) h ;u -u h ,v- v h , v) 

< A(u -u h ,v- v h , v) + j(v h ) - j(v) - B(u h , v h , <p h ;u -u h ,v- v h , v) 

< A(u -Uh~ (uh - u h ),v - v h - (v h - v h ),u- u h ) +j(v h ) ~ j(v) 

- B(u h ,Vh,(f>h;u ~ Uh ~ {Uh ~ u h ),v -v h ~ {Vh ~ v h ),v- u h ). 

Note that the variational inequalities are identities when restricted to the 
^-variable. The claim follows by setting eh = Uh — Uh and e~h = Vh — Vh- □ 

Simplifying the right hand side along the lines of [2j leads to a gradient 
recovery scheme in the interior with a residual type estimator on the bound- 
ary. With a straight forward modification of [8] , also a method purely based 
on residual type estimates could be justified. 

For 1 < p < oo and < 8 < 1, define 

G p ,s(x,y) = \y\ 2 u{x,yf- 2 = \y\ 2 [{\x\ + \y\f{l + \x\ + lyl) 1 -^ 2 

whenever \x\ + \y\ > and otherwise. As in [2], our analysis will be based 
on the following consequences of the monotony and convexity properties of 
G Pj s- 
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Lemma 6. Assume that f2 is connected. Let q be a continuous linear form 
on W 1,P (Q) with M D kerg = {0}, where R is identified with the space of 
constant functions on £1. Then for any 1 < p < oo there exists Cp = 
Cp(p, q, Q) > such that for all a > and u G W 1,p (fi), 

/ G p<s (a,u) < C P ( G p> s(a,q(u)) + / G Pj(5 (a, |Vn|) 
in V in 

Proof. Cf. [2], Lemma 4.1 and its generalization in Remark 4.3. □ 

Lemma 7. For any d, k G N i/iere zs Cs = Cs(p, d, A;) > suc/i i/iaf /or a// 
01,02, • • • ,Ofc G R d 

k 3-1 fe-1 
y^y^ G p^( a i' a i ~ a «) ~ V] min G p S (a m , a j+ i - aj). 

1 — ' ' — ' 1 — ' Km<k 
j=l 1=1 j=l ~ ~ 

Proof. Cf. [2], Lemma 4.2 and its generalization in Remark 4.3. □ 

Even though Lemma and Lemma [9] hold for any 1 < p < oo with minor 
modifications of the proofs (see [2] for a similar discussion), we will from 
now on concentrate on the range 2 < p < oo relevant to our transmission 
problem. 

Definition 3. Let z G D be a node of the triangulation Th and tp z G W^ P (Q) 
the associated nodal basis function. Let lo z = {x G Q : ^p z (x) > 0} be the 
interior of the support of ip z . The interpolation operator it : W l,p (Q) — > 
Wl' p {VL) is defined as 



TTU = 

z&D 

Lemma 8. Let £ ? = {I G £h : I = K{ D JKj /or some Ki,Kj C Given 
Uh G W^' p (r2), Zei [<9^ £ n/j] ; denote the jump of the normal derivative across 
the inner edge I of the triangulation. Then, if v G W 1,p (fJ) and G %h, the 
following estimate holds: 

G Pj s(Vu h ,h] ( 1 (v - nv)) + / G PjS (Vu h ,V(v -ttv)) 
k Jk 



/ Gp s^Vuh, Vv) + V] min / G p ,s(Vu h \ K ', [9 V£ u h ],) \ 
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Proof. The proof is a modification of [2], Lemma 4.3. Concerning the first 
term on the left hand side, the convexity of G p ^ in its second argument (a 
"triangle inequality") and enlarging the domain of integration leads to 



/ Gpj(Vuh,hjJ-(v -irv)) = / G PtS (Vu h , h K( v - v z)Vz) 
Jk Jk zeDnR 

~ / G P)5 {Vu h ,h^(v - v z )tp g ) 

zDnR 

^2 / Gl PtS iyu h \ K ,K^{v -v z )(p z ). 



zeDnK 

< 



zeDnK 



As G P! s(Vuh\K,') is increasing and \ip z \ < 1, Lemma[6]with q(u) = J (p z u 
implies 



Gp^Vut^Kih^iy - v z )<p z ) < / G p>s (Vu h \ K , hx(v - v z )) 

(3) < C P I G p>s (Vu h \ K ,V(v-v z )) 

= C P I G p , s (Vu h \ K ,Vv) 

for every term in the sum over z 6 D n K. To replace the constant Vuh\x 
by Vuh, we repeatedly apply the usual triangle inequality and the convexity 
of G Pt s to obtain 

G Pj s(yu h \ K ,Vv) 

< G PiS (Vu h \ K , \Vv\ + \Vu h \ K - Vu h \) 

= (\Vv\ + \Vu h \ K - Vu h \f(\Vu h \ K \ + \Vv\ + \Vu h \ K - Vu h \ f p -V 
x (1 + \Vu h \ K \ + \Vv\ + \Vu h \ K - VuhD^-W-Q 

< (\Vv\ + \Vu h \ K - Vu h \) 2 (\Vv\ + 2(|V« fc | + IVufcl* - Vu h \)) 5 ^ 
x (1 + |V«| + 2(|V« h | + \Vu h \ K - Vufcl))^^ 

< G PtS (X7u h , \Vv\ + \Vu h \ K - Vu h \) 

< G PjS (Vu h , Vv) + G PiS (Vu h , Vu h \ K - Vu ft ). 

Altogether 

/ G Pi s(S! f Uh\K,h~K(v-'Kv))< ^2 / {G P ,s(Vu h , Vv) + G p ,s(Vu h , Vu h \ K - Vu h )} . 
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Let u z = K\ U • • • U Kk- Applying Lemma [7] with aj = Vu^Kj, 1 < J < A;, 
leads to the asserted bound for the first term. For the proof, note that the 
conormal derivatives of the piecewise linear function Uh are determined by its 
boundary values on the corresponding edge. But Uh G W^ P (Q) C W 1,p (£l), 
so the restrictions from both sides have to coincide, and the conormal deriva- 
tive does not jump: aj — a,j-\ = [dv e u h\Kjr\Kj-i\- 
As for the second term, let c = t^t J k v. Because 

/ G p . s (Vu h ,V(v-irv))< [ G P 4X7u h ,Vv)+ [ G p , s (Vu h , V(nv - c)) 

JK JK JK 

by convexity and the triangle inequality, it only remains to consider the 
second term J K G Pt s(Vuh, V(nv — c)). The inverse estimate 

\V(irv - c)| < jf h^lnv - c\ 

for the affine function irv — c and Jensen's inequality show 

/ G p . s (Vu h ,V(7TV-c)) < [ [ G PtS {Vu hl h- K l {-Kv-c)) 

JK JK l A I JK 



G PtS (Vu h , h K l (irv - c)). 

K 



However, as before 



/ Gp^Vu^hjl^v-c)) < / G Pi s{Vu h ,h K 1 (v-7Tv))+ G Pt s(yu h ,hjl(v- 

JK JK JK 

and the first term has been considered in the first step of the proof. Lemma [6] 
with q(u) = J K u also bounds the final term by f K Gp^Vu/j, Vv). □ 

Lemma 9. For any e > 0, u h e W^ P (Q), v G W l > p {£l) and f G L P '{Q), 
I f(v-n V ) < Ce [ G P)& {Vu h yv) 

JQ JQ 

zeDKCZJz Jk 

Here, fx = jj^ J K f ■ If f G W 1,p (O), the second term may be replaced by 
C ( £ )E E / G plyl {\Vu h \ p -\h 2 K Vf). 
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Proof. We adapt the proof of [2], Lemma 4.4. Let K C u z such that 
l^ u /il_ft-| = max K'cuj z Applying the inequality from Remark [2k) 

for some e > and C(e) = Cp max{e _1 , e 1 /( 1 -p)]. ) 

f(v-irv) = ^2^2 / h K{f ~ fK)h^{v -v z )ip z 

< C p l C{e)Y,Y. /(|V^|^r 1 + ^|/-/K|) p '- 2 ^|/ 

z&DKCTT K 



+ £ ^2 / d Vn ^l/rl + ~ v z \(p z ) p 2 h K 2 \v - v z \ 2 tf 
< Cp l C{e)Y,Y. [ G P',i(\Vu h \ & \ p -\h K (f-fK)) 

zeDKCZJ z ^ K 
zeDKCTT JK 



[V ~ V z )if z ), 



because YIkcjd Ik f K ( v ~ v z) L Pz = 0. However, by our choice of K and 
because p' < 2, 

[ G pf>1 (\Vu h \ K \ p -\hK(f-fK))< [ G p , A (\Vu h r\h K (f-f K )). 

J K JK 
If / G W^P'ifl), Lemma [6] with q(u) = J R u gives: 

/ GV,i(!V^r\M/ - f K ))<C P [ G p , A (\Vu h r\h K \7f). 

JK JK 

Concerning the Gp^-term, equation ([3]) in the proof of Lemma [8] shows that 
it is dominated by e G Pi s(yuh\ K , Vi>), which in turn was bounded by 

£ i G pS (Vu h ,Vv) +e y _min / G p ^(Vm^, [d Ve u h ]X 

□ 

In order to define the a posteriori estimator, we still need to introduce 
some notation. For any z € D, denote by Kj z G Th, 1 < j < N z , the 
triangles neighboring z in the sense that co z = Uj=i Kj,z- To each Kj z we 
associate a weight factor aj jZ > normalized to X/i=i a j,z = 1- 
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Definition 4. Given Uh G W^ P {£1), define the gradient recovery 

N z 

G h uh = ^2(GhVh)(z) <Pz, (G h Vh)(z) = ^2a jtZ Vu h \ Kj!Z - 
z£D j=\ 

The following theorem states our reliable, but presumably not efficient 
a posteriori estimate. 

Theorem 3. Let f G L' p ' (Q) and denote by (e, e, e) the error between the 
Galerkin solution (u,h,Vh,(ph) £ Y p anc ^ ^ e ^ rue so ^^ on {u,v,(f)) G Y p . If 
T s ^ 0, assume that Vu|r s G L P (T S ). Then 



II?,. < 



where 

'Igr 



/ G P,s(^u h ,X7u h - G h u h ), 
r l)= E f K G P >,i(\Vuh\ p -\h K (f-f K )), 



KeT h ' 



4=Y, 1 W d s{ v <t>h + (1 - K)(« A |en + «fc - «o)}llj 



Il2(o 



^ = E l \\e(Vuh) dvUhlr.Wh® + lbll^-i, 2(rs) 

Iff e W l *\tt), we may replace n) by ^KeT h Ik G p' ' ,i(\™h\ p -\ h 2 K V f). 
Proof. From Lemma [5] we know that for all (e^, v^) G Y? 

II ~ IIP I |2 i II I i — 1 1 2 , II i|2 

ii ii yp ~ i i(i,u,pj n 'V^' 2 (an) 11 'V _ 2-- 2 (aa) 

^ / /( e - e h) ~ / ^( V ^) d uUh\dK (e - e A ) 

+ / 3(14 + Vh\ ~ \v\) ~ {v ~ v h ,V4>h + (1 -£){uh\dsi + Vh -uo)) 
+ (*o - W{uh\dn + v h - u ) - (K! - l)4> h , (e - e h )\ dn + e - e ft ), 
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with 2v = e + V _1 (l — /C)(e|an + e). The first two terms are mainly going 
to give the gradient recovery in the interior, the fourth term the error rjs of 
constructing the Steklov-Poincare operator, while the remaining terms add 

Up tO 7]Q. 

Concerning the first term: 



[ f(e-e h ) < I G PtS (Vu h ,Ve) 

Jn KeTh Jk 

+C{e)Y, [ G p/A (\Vu h r\h z (f-f z )) 
+eVV_min / G pS (Vuh\K',[ d i's u h]i) 

£ £ \ e \ 2 (l,u, P ) + C ( £ ) 7 l} + £ YjYj / G Pj(™h\K>, [dueUhW 

z eDie£* Knl ^ Juj * 



GhUh is continuous across any interior edge I, so that [d u Uh]i = [dvUh—GhUh]i 
and 

_min / G Pt s(Vu h \ K >,[d u u h - G h u h }i) < G Pt s(Vu h ,Vu h - G h u h ). 
Therefore, 

/ f(e-e h ) < e\e\ 2 {lj ~ :p) + C(e)r? f + e ^ ^ / G p , & {Vu h ,[d u u h - G h u h ]i) 
~ £ \ e \ 2 (i,u, P ) + c ( £ Wf + e V / G PtS (Vu h ,Vu h - G h u h ) 

= £ \ e \tl,u, P ) + C ( £ H + £ Vgv 

Concerning the second term, let 

A t = Q(Vu h \ KlA ) d u u h \ Kl l - g(Vu h \ Kli2 ) d u u h \ Kl/2 , 

where again I C Ki\ n Kip, and the unit normal v points outward of K[ \. 
Therefore 

-(G'u h ,e-ire) = - ^ / q(Vu h ) d v u h \ dK {e - ire) 

Ker h JdK 

^2 / A t (e - ire) - ^ / g(Vu h ) d u ii h \i (e - ire). 



i<£dn " 1 icdn • 



19 



Repeating the analysis of [2], Theorem 5.1, with the help of Lemma [8] gives 
J Me - vre) < r, 2 gr + e(\e\\ xM + r] 2 gr ). 



lt£dQ' 



Thus 



IP <? I |2 i II I i ~||2 , n ||2 

~ \ e \(iuv) + \\ e \dn+e\\ i„ + e i ., 

~ *7/ + e ( 7 5r + ! e l(l,«,p)) + ^sr + e (M(l,u h ,p) + 

+ / {-^(Vn^) ^-Uh|r s (eh - e) + g(\e h + v h \ - \e + v h \)} 
Jt„ 

£>(Vu/i) ^%|an ((e - 7re)| 9 n + e - e h ) 

+{t - W(u h \ d n + Vh~ uo) ~ (£■' ~ (e - ne)\ dn + e - e h ) 
-(u- v h ,V4>h + (1 - + Vh ~ uo)}- 

We bound the second, third + fourth as well as the final line individually. 
Cauchy-Schwarz and Young's inequality allow to estimate the last term by 

and the latter by rjg (cf. [3.]). The third and fourth lines are estimated by 

(*. 131) 

||-£>(V%) 9 v %+i o -W(%|an-H)/ l - , uo)-(^ / -l)0h|| w ,-^a (an) II ( e - 7re )!an+e|| H/ i 

which lead to 779, where we have choosen = 0, i.e. Vh = Vh- Finally, using 
the triangle inequality, the second line is simplified as follows: 



{-g(Vu h ) d v u h \r s {e h - e) +g(\e h + v h \ - \e + v h \)} 
< / {-Q(Vu h ) d v u h \r a (e h -e)+g\e h -e\} 
{?(V%) d u u h \ Ts e + g\e\} 



< WoCVuh) d v Uh\rA\ 1 9, Jlell 1 , + llgll 1 \\e\\ 1 9 , 

We may use the Cauchy-Schwartz inequality and the inverse inequality, lead- 
ing to ri g . 

□ 
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7 Numerical results 

~I 2 

With the subset A h of W£' (F s ) given by 

Aft = Wh 6 ^|' 2 (r s ) : |<7 fc (x)| < 1 a.e. on T s }, 

we can define an Uzawa algorithm for solving the variational inequality 

analogously to [9]. In order to introduce this algorithm, let Pa be the 

— i 2 

projection of W£' (T s ) onto A/j, i.e. for every nodal point of the mesh 
T h \ Ts holds 5 P A (5) = sup{-l, inf(l, 5)}. 

Algorithm 1 (Uzawa). 
1. Choose o\ G A ft . 

g. For n = 0, 1, 2, . . . /md ug) G X£ suc/i that 

(G'v%,u h ) + (S h (v%\dn + v%),u h \ dn + v h ) + / ga%v h ds = \ h (u h ,v h ) 

Jr s 

for all (u h ,v h ) G X v h . 
3. Set 

a n h +1 = P A K + P gv n h ), 
where p > is a sufficiently small parameter that will be specified later. 

4- Repeat with 2. until a convergence criterion is satisfied. 

In our first example the model problem is defined on the L-shape with 
n = [~h l] 2 \[°> il 2 ' QC = ffi2 \^- The friction part of the interface is T s = 
H,-l)(l,-l) U (-|,-|)(-|,|), see Figure HJ 

In this example we choose g(t) = (e + t) p ~ 2 , with p = 3 and e = 0.00001. 
Our volume and boundary data are given by / = and no = r 2//3 sin \ (}p— § ), 
io = 9 u uo\qq. The friction parameter is g = 0.5, leading to slip conditions 
on the interface. We have applied the Uzawa algorithm as introduced above 
with the damping parameter p = 25 to solve the variational inequality. 
The nonlinear variational problem in the Uzawa algorithm is then solved by 
Newton's method in every Uzawa-iteration step. 

In Table Q] we give the degrees of freedom, the value Jh(uh,Vh) and the 
error measured with the help of J, i.e. 5 J = Jh(uh,Vh) — J(u,v), where we 
have obtained the value J(u,v) by extrapolation of Jh(uh,Vh)- Due to the 
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slip condition, we need only a few Uzawa steps. But as a consequence of 
the degeneration of the system matrix, due to the nonlinearity, the iteration 
numbers for the MINRES solver, applied to the linearized system, are very 
high, leading to large computation times. The convergence rate aj is sub- 
optimal, due to the presence of singularities, in the boundary data as well, 
as due to the change of boundary conditions. 



n 



Figure 1: Geometry and interface of the model problem 
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Table 1: Convergence rates and Uzawa steps for uniform meshes (Example 
1) 

In our second example we use the same model geometry as before (see 
Fig. H|). Here we choose the friction boundary T s = 0. Therefore our model 
problem reduces to a non-linear p-Laplacian FEM-BEM coupling problem, 
where we can prescribe the solution. 

In this example we choose g(t) = (e + t) p ~ 2 , with p = 3 and s = 0.00001. 
We prescribe the solution by U\ = r 2 ' 3 sin|(<£ — ^) and 112 = 0. Then 
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the boundary data uo,to and volume data / are given by uq = Ui|r> to = 
g{\Vui\)d u ui and / = — div(fj(|Vui|)Vui). 

In the following we give errors in the || • ||vm>p(C) norm and in the quasi- 
norm \u - u h \ Q = \\u - u h \\(i )UhjP y 

In Tab. we give the errors, convergence rates, number of Newton it- 
erations It N ew ton and the computing time for the uniform h- version with 
rectangles. We observe that the convergence rate in the quasi- norm \ ■ \q 
is better than in the || • ||^i,3(Q)-norm. The number of Newton iterations 
appears to be bounded. 

In Tab.[3]for the uniform h- version with triangles, we give the errors, con- 
vergence rates, error estimator 77, efficiency indices 5 u /r] for the || • ||wi.3(n) _ 
norm and 5 q /r] for the | • |g-norm, number of Newton iterations and the 
computing time. Again, here we observe that the convergence rate in the 
quasi-norm | • \q is better than in the || • ||yj/i,3(£})-norm and the number 
of Newton iterations is bounded. The efficiency index b u jr\ appears to be 
constant, whereas the efficiency index 5 q /r] appears to be decreasing. 

Tab. H] gives the corresponding numbers for the adaptive version, using a 
blue-green refining strategy for triangles and refining the 10% elements with 
the largest indicators. Here we observe that the convergence rates for both 
norms are very similar and that both efficiency indices are bounded. 



Figure [2] give the errors for all methods in the 



•norm and 



the I • |q quasi-norm together with the error indicators for the uniform and 
adaptive methods. 

Figure [3] presents the sequence of meshes generated by the adaptive 
refinement strategy. We clearly observe the refinement towards the reentrant 
corner with the singularity of the solution. 



conf-uni-h-3 ------- 

adap-residual-h-3, thela=0.90 

conl-uni-h-3-indic □ 
adap-residual-h-3, thet a =0.90- indlc ■ 






conf-uni-h-3 - 
adap-residual-h-3, iheta=0.90 

conl-uni-h-3-indic □ 
adap-residual-h-3, fheta=0.90-indic 
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Table 2: Errors, convergence rates (Example 2, uniform mesh with rectan- 
gles) 
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23 


144.2 


12545 


0.0611994 


-0.168 


0.0283778 


-0.246 


0.236 


0.260 


0.120 


24 


608.7 


49665 


0.0486160 


-0.167 


0.0203130 


-0.243 


0.184 


0.265 


0.111 


24 


2530. 


197633 


0.0386151 


-0.167 


0.0145686 


-0.241 


0.144 


0.269 


0.102 


24 


11000 



Table 3: Errors, onvergence rates, estimator r/, reliability S u /rj and S q /rj 
(Example 2, uniform mesh with triangles) 
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DOF 


\\ u ~ u h\\l,3 


a 


\u ~ u h \q 




Qt 




n 


5n If) 


Sq/r] It New 


t(s) 


21 


0. 


.1945908 






0. 


.1510064 






1 


,027 


0.190 


0, 


147 


22 


0.196 


32 


0. 


.1602214 


-0. 


.461 


0. 


.1205155 


-0. 


535 


0. 


.804 


0.199 


0. 


150 


22 


0.332 


54 


0. 


.1275298 


-0. 


.436 


0. 


,0918131 


-0. 


.520 


0. 


,603 


0.212 


0. 


152 


22 


0.648 


93 


0. 


.1019990 


-0. 


.411 


0. 


.0699054 


-0. 


.501 





,442 


0.231 


0. 


.158 


22 


1.132 


152 


0. 


.0821754 


-0. 


.440 


0. 


,0540462 


-0 


.524 


0. 


.325 


0.253 


0, 


166 


23 


2.000 


249 


0. 


.0679251 


-o. 


.386 


0. 


.0449420 


-0 


.374 


0. 


.246 


0.276 


0, 


183 


23 


3.352 


400 


0. 


.0558447 


-0. 


.413 


0. 


.0369614 


-0. 


.412 


0. 


,190 


0.294 


0. 


.194 


23 


5.700 


625 


0. 


.0439784 


-0. 


.535 


0. 


.0277857 


-0 


639 


0. 


.148 


0.297 


0. 


.188 


24 


9.896 


986 


0. 


.0352491 


-o. 


.485 


0. 


.0217361 


-0. 


539 





.116 


0.305 


0. 


188 


24 


17.45 


1528 


0. 


.0279287 


-o. 


.531 


0. 


.0167409 


-0. 


.596 


0. 


.091 


0.308 


0. 


.184 


25 


31.16 


2322 


0. 


.0222760 


-0. 


.540 


0. 


,0129489 


-0. 


.614 


0. 


.071 


0.312 


0. 


181 


25 


53.98 


3620 


0. 


.0177640 


-0. 


.510 


0. 


,0102552 


-0. 


.525 


0. 


,056 


0.316 


0. 


182 


25 


106.7 


5544 


0. 


.0142059 


-o. 


.524 


0. 


.0080233 


-0 


.576 





.044 


0.320 


0. 


181 


25 


205.3 


8449 


0. 


.0112965 


-o. 


.544 


0. 


.0063426 


-0 


,558 


0. 


.035 


0.322 


0. 


181 


26 


422.4 


12810 


0. 


.0090396 


-o. 


.536 


0. 


.0050706 


-0 


.538 





.028 


0.325 


0. 


183 


26 


1060. 


19222 


0. 


.0072288 


-0. 


.551 





.0040370 


-0 


.562 


0. 


,022 


0.329 


0. 


.184 


26 


2400. 


29006 


0. 


.0057984 


-0. 


.536 


0. 


.0032478 


-0. 


529 


0. 


,018 


0.333 


0. 


186 


27 


5460. 


43593 


0. 


0046615 


-o. 


.536 


0. 


.0026230 


-0 


,524 





.014 


0.337 


0. 


190 


27 


13000 



Table 4: p-Laplacian (adaptive), convergence rates, estimator 77, reliability 
5 u /r] and 5 q /r) 
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Figure 3: The first 12 meshes generated by the adaptive refinement alg 
rithm 
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